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Abstract 

We introduce a model described in terms of a scalar velocity field on a Id lattice, 
evolving through collisions that conserve momentum but do not conserve energy. 
Such a system posseses some of the main ingredients of fluidized granular media and 
naturally models them. We deduce non-linear fluctuating hydrodynamics equations 
for the macroscopic velocity and temperature fields, which replicate the hydrody¬ 
namics of shear modes in a granular fluid. Moreover, this Landau-like fluctuating 
hydrodynamics predicts an essential part of the peculiar behaviour of granular flu¬ 
ids, like the instability of homogeneous cooling state at large size or inelasticity. 
We compute also the exact shape of long range spatial correlations which, even 
far from the instability, have the physical consequence of noticeably modifying the 
cooling rate. This effect, which stems from momentum conservation, has not been 
previously reported in the realm of granular fluids. 


1 Introduction 

Since the seminal paper of Einstein |l|, it is well known that the fluctuating behaviour 
of systems at the mesoscopic level reflects the hectic microscopic dynamics beneath. 
While the equilibrium behavior of mesoscopic fluctuations has been investigated and 
understood in detail 12,3], a big effort is still being carried out to explore the fluctuating 
properties of non-equilibrium media |4|. These are known to lead, in great generality, to 
the emergence of spatial correlations and pattern self-organization [5j[6]. In this context, 
the crucial task of connecting microscopic and mesoscopic dynamics is considerably 
simplified when there exists a separation of scales, which makes it possible to introduce 
slow fields evolving under the so-called hydrodynamic equations [7j. 
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An important class of systems exhibiting patterns includes two types of complex 
fluids: active matter j8,9|, such as bacteria or birds, and fluidized granular materials 110 . 
Interestingly, active and granular matter are often associated 111 14]. They are not only 
relevant for applied and biomedical sciences, but also offer fascinating challenges for 
kinetic theory 15 . Indeed, the lack of energy conservation in the microscopic dynamics 

In granular and active fluids, 
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makes them intrinsically out-of-equilibrium systems 
the spectacular emergence of spatial patterns, particularly in vectorial fields such as 


momentum or orientation, is often understood in terms of hydrodynamic equations 17 


Furthermore, a relevant role is played by fluctuations, as an inevitable consequence of 
the relative small number of their elementary constituents 
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One of the most intensively investigated states in the realm of granular fluids is the 

Therein, the granular temperature decays in 
whereas the system remains spatially homogeneous. Re- 


homogeneous cooling state (HCS) 119| 20 
time following Haff’s law 


markably, the HCS is the reference state for the hydrodynamic description of granular 
fluids but is unstable: for large enough inelasticity or system size, the scaled fluctuations 
of the transverse velocity increase (shear instability) and eventually density inhomo¬ 
geneities arise (clustering instability) |: 


. This instability for large system sizes makes 
it relevant to look into the finite size corrections to the physical quantities, like the cool¬ 
ing rate, since a “thermodynamic limit” in which the system size is infinitely large cannot 
be taken without simultaneously making the inelasticity vanish. Notwithstanding, and 
to the best of our knowledge, these finite size corrections have only been investigated for 
a system of smooth inelastic hard spheres very close to the shear instability 
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In deriving mesoscopic transport equations from microscopic rules, analytical results 
are needed and simple models are good candidates for this (24|25|. In this paper, we study 
a Id lattice model which implements two main ingredients of granular fluids: inelastic 
collisions and momentum conservation. Given the simplicity and appealing physical 
picture of the model, this novel approach may help to improve our current understanding 
of the complex behaviour of granular fluids. To start with, we recover some of the 
main features of granular fluids: more specifically, the large size and inelasticity shear 
instability. In addition, we are able to compute exactly the shape of velocity correlations, 
which allows us to, first, extend the known results of fluctuating hydrodynamics 26 


and, second, obtain the finite size corrections to the cooling rate. The latter has a non¬ 
trivial dependence on the inelasticity and system size: for a given inelasticity, it changes 
sign at a certain value of the system size that is smaller than the one corresponding to 
the shear instability. 

Finally, we would like to stress that momentum conservation is a physical constraint 
that certainly has a recognised role in the appearance of long range spatial correlations |6, 
27 . However, since it complicates the description and the derivation of exact results, 
it is rarely considered in its entirety. Here, starting from the microscopic dynamics, 
we are able to rigorously derive the mesoscopic equations that describe the average 
and fluctuating behaviour at the hydrodinamic scale, taking into account momentum 
conservation in full. 
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2 Microscopic equations of the model. 


Fluctuating hydrodinamics in linear and nonlinear lattice diffusive models have been 
extensively studied in recent years, both in the conservative 28 -31 and in the dissipative 
cases |32j-[34|. Inspired by 35 , we consider a Id lattice with N sites and given boundary 


conditions, either periodic or thermostatted, depending on the situation of interest. At 
a given time p, each site l possesses a velocity vi tP and the total energy of the system is 


E p = vf p . In an elementary time step of the dynamics, with a probability discussed 
below, a pair of nearest neighbors (l, l + 1) collides inelastically and evolves following the 
rule (0 < a < 1) 


^i,p +i — ^ l,p (1 3“ &)Ai p /2 , ^i+i,p+i — ^(+i,p “I - (1 “I - oc)Ai >p / 2, (1) 

having defined the relative velocity 

A/ : p — vi,p ^(+i,p- (2) 

Momentum is always conserved, vi tP + f;+i,p = v l,p +1 + vj+i,p+i, while energy, if a / 1, 
is not: v 2 p+1 + v? +lp+1 - vl p - vf +l p = (a 2 - l)A 2 p /2 < 0. 

The definition of the model implies that there is no mass transport, particles are 
at fixed positions and they only exchange momentum and kinetic energy. We are also 
disregarding the so-called kinematic constraint in [35] , namely a colliding pair is chosen 
independently of the sign of its relative velocity. This can be understood as the velocity 
of the particles representing not their motion along the lattice axis but rather along a 
transversal one: in fact, the hydrodynamics derived here replicates transport equations 
for granular gases in d > 1 restricted to the shear (transverse) velocity held, see Fig. [I] 
In the context of granular fluids, the model may be physically motivated as follows. 
We start from a d > 1 system that has been divided into “slabs” that are perpendicular 
to the lattice direction. Specifically, each particle on the lattice represents one slab. 
In this sense, the parameter a that appears in the collision rule @ should not be 
confused with the usual restitution coefficient defined in granular media, since here a 
stands for an effective inelasticity for the collisions between slabs. The connection with 
a “real” granular fluid should be done at the level of the cooling rate that appears in 
the hydrodynamic equations, see Section [3] 

Now we write down the evolution equation for the velocities. At time p, the probabil¬ 
ity that the nearest-neighbours pair at sites (l, Z+l) collide is assumed to be P^ p oc | A/ )P | /3 . 
Then, 

vi,p+ 1 Wp = Jl,p "F 3i,p = (1 d o / -)Ai p dy p ^i/2, (3) 

which is a discrete continuity equation for the (conserved) momentum, with ji p being 
the momentum current, that is, the flux of momentum from site l to site l + 1 at the 
time step p. Therein, 6 yp j is Kronecker’s 5 and y p is a random integer which selects the 
colliding pair with probability Pi )P - The evolution equation for the energy is obtained 
by squaring 0, 

v i,p+ l ~ v i,p = ~Ji,p + Ji-i,p + di,p- (4) 
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t'i,p+i - *i,p (1 t ®)^i,p/2 

t'i+i,p+i = %i,p+(l+o)A/.p/2 



Figure 1: The model is defined on a lattice, each site being characterized by a velocity 
vi and standing for a fluid slab. The dynamics proceeds via nearest-neighbour collisions 
as defined in the text, in which part of the energy is dissipated. This model resembles 
a sheared granular system at the mesoscopic level. To be specific, we show the sketch 
that corresponds to periodic boundary conditions. 
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Again, we have defined an energy current J^ p from site l to site l + 1 and the energy 
dissipation d^ p at site l as 

a 2 — 1 

Ji,p = ( v i,p + v i+i,p)ji,pi di, P = ^ (dy pt iAi p + 6y py i-iAi_ l p ) < 0 . ( 5 ) 

The sink term only vanishes in the elastic case a. = 1. 

Stochastic simulations of the model are done as follows. At any Montecarlo step 
p , one site l is picked with probability Pi )P oc |A/ p |^ and particles l and l + 1 collide 
following the microscopic rules ([Tj) . The simplest choice for Pi iP corresponds to /3 = 0, 
all the pairs are chosen with uniform probability, Pi tP = 1/L, in which L is the number 
of pairs. This is often called in the literature the model of inelastic Maxwell molecules 
(MM) (35 1 . Note that L is basically equal to N but depends on the boundary conditions: 
for periodic boundary conditions, it is L = N, but if we consider the system coupled to 
two extra sites 0 and N + 1, which introduce the appropriate boundary conditions, it is 
L = N + 1. The periodic boundary conditions, sketched in figure [lj correspond to the 
free (undriven) evolution of the system and if l = N it is the pair ( N , 1) that collides. 

In the following, we discuss the hydrodynamic limit, fluctuations and correlations 
for a particular choice of Pi tP . Specifically, we consider the case of MM: such a choice is 
dictated by the will of simplifying the presentation and making clear the essential points. 
We postpone a more complete and general discussion to a more technical and detailed 
paper, in preparation. The theoretical results are compared to the numerical simulations 
described above. A large enough value of L, which is indicated in the figures, has been 
used to ensure the hydrodynamic limit, and we have averaged over 10' 5 realizations of 
the stochastic dynamics. Aside from MM, results for HS (/3 = 1) are also shown in a 
few, clearly marked, cases. 

3 Hydrodynamic limit: average equations and fluctuations 

Let us define, as usual, the following local averages over initial conditions and noise 
realizations: = ( vi jP ), E) )P = {vf p } and T[ p = E^ p — u 2 p . Their evolution is obtained 

under a series of assumptions. With the choice of MM, y p is an uniform distributed 
random integer, namely (5y p ,z) = 1 /L, with L being the number of nearest neighbor 
pairs. In addition, when considering the average dissipation at site l , there appears 
moments like (vi )P vi~_ ti, p ). To the lowest order, we assume that neighbouring velocities 
are uncorrelated, that is, (vi )P vi±\ tP ) = ui tP ui±\ tP (see Appendix A). 

Now we assume that ui }P and E^ p are smooth functions of space and time and 
introduce a continuum (“hydrodynamic”) limit (HL) by defining macroscopic scales: 
Ax ~ L ~ 1 and At ~ L~ 3 . Each spatial derivative introduces thus a factor L~ l in 
the continuum limit: therefore, the difference between the current terms in the balance 
equations is of the order of L~ 3 . On the other hand, the dissipation goes as (1 — a 2 )^ 1 , 
which makes it useful to define the cooling rate as (see Appendix |Aj) 

v = (1 — a 2 )L 2 . (6) 
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It is natural, on the scales defined by the HL, to define the mesoscopic fields u(x,t), 
E(x, t ) and T(x, t), as well as the average mesoscopic currents j(x, t ) = lirri/^oo L 2 (jy p ) = 
—d x u(x,t), J(x,t) = lim^oo L 2 ( J; iP ) = —d x [u 2 (x,t) + T(x,t)] and the average meso¬ 
scopic dissipation of energy d(x,t ) = lim^o Q L 3 (dy p ) = —uT(x,t), which depends only 
on the granular temperature but not on the average local velocity u(x,t), as expected 
on a physical basis. After computations detailed in Appendix [Aj we get the HL of (]3]) 
and Q, which are dtu(x,t ) = —d x j(x,t) and dtE(x,t ) = —d x J(x,t) — uT(x,t) respec¬ 
tively. Then, we can write the average hydrodynamic evolution equations 

d t u = d xx u , d t T = -uT + d xx T + 2 ( d x u ) 2 . (7) 


over the length and time scales defined above. Note that, here, we have substituted 
1 + a by 2, because a 2 = 1 — vL~ 2 , and we have already neglected L^ 1 terms. These 
equations must be supplemented by appropriate boundary conditions for the situation 
of interest. Technical details are deferred to a later paper. 

Let us consider the fluctuations of the microscopic currents and dissipation, that 
is, ji, p = ji, P + €i, p , Ji, p = Ji, p + Vl,p, an d di tP = di tP + 9i iP . Tilde variables correspond 
to a partial average: they are averaged over the fast variables y^ p but not over the 
slow ones vi )P . Thus, for example, ji tP = (1 + a)Ai tP /2L. It is clear that this choice 
guarantees that all noises £/ iP , r]l,p and Q^ p have zero average. The noise correlations read 
(££') ~ A^S(x — x')5(t — t'), (yy'} ~ A v 5(x — x')5(t — t') and (99') ~ Ag5(x — x')5(t — t') 
with amplitudes A g = 2L~ 1 T(x, t), A p = 4 L~ 1 T(x,t)[T(x,t) + 2 u 2 (x,t)\ and Ag = 
3L~ 3 v 2 T 2 (x,t) (see Appendix [b]). In the above relations, we have used the notation 
£ = £(x,t) and £ = x',t and similar notations for y,y',9,9'. Thus, the currents 
noises are delta correlated in space and time, and their amplitudes scale as L~ l with 
the system size L. On the other hand, the noise of the dissipation is subdominant with 
respect to the moment and energy currents, its amplitude being proportional to L -3 , and 
therefore it is usually negligible. Gaussianity for these noises can be easily demostrated, 
see 32 . Interestingly, being in the presence of two fluctuating fields, correlations between 
different noises appear. Theoretical predictions for noise correlations, amplitudes and 
Gaussianity have been succesfully tested in both MM and HS simulations, see Appendices 

HandO ' 


4 Solutions, HCS and instabilities 

Here we focus our attention on the case of spatial periodic boundary conditions and an 
initial “thermal condition”: vy o is a random Gaussian variable with zero average and 
unit variance, that is, Tj q = T(x, 0) = 1. Starting from this condition, the system is 
expected to typically fall into the so-called Homogeneus Cooling State (HCS), in which 
the velocity and temperature fields remain spatially uniform, and the temperature decays 
in time. Indeed, the solution of the average hydrodynamic equations 0 reads 

u(x, t) = 0, T HC s(x , t) = T(t = 0)e~ ut . (8) 


6 



The exponential decrease of the granular temperature is typical of MM, where the col¬ 
lision frequency is velocity-independent. It replaces the so-called Haff’s law which was 

because T oc T 3 / 2 


originally derived in the HS case, where Thcs ~ t 


-2 
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The HCS is known to be unstable: it breaks down in too large or too inelastic sys¬ 


tems 36 . In our model and in the hydrodynamic limit, this condition is expected to be 
replaced by a condition of large u. The stability is studied by introducing rescaled fields 
U(x,t) = u(x, t)/yfTncs{t) and T = T(x,t)/THcs(t) and by linearizing the hydrodi- 
namic equations near the HCS, i.e. T(x,t ) = Tncsif) + ST(x,t) and U(x,t) = 5U(x,t). 
The analysis of linear equations becomes straightforward by space-Fourier-transforming, 
which gives 


dtSU(k,t ) = 


v — 2k 2 


SU(k, t ), d t 8T(k, t ) = —k 2 5T(k, t). 


(9) 


Therefore, 5U is unstable for wave numbers that verify v — 2 k 2 > 0. In the continuous 
variables we are using, the system size is 1, so that the minimum available wavenumber 
is k m i n = 27T. Thus, there is no unstable mode for v (lengths) below a certain threshold 
v c (L c ), with 

v c = 8ir 2 , L c = 2nV2(l- a 2 y 1/2 . (10) 

On the contrary, for u > v c (L > L c ), the HCS is unstable and modes with wave 
numbers verifying k < \Jv/2 increase with time. This instability mechanism is identical 
to the one found in granular gases for shear modes |26j. Theoretical predictions and 
simulations results perfectly agree, as plotted in Fig. [2] It is important to stress that 
the amplification appears in the rescaled velocity U(x,t ) and not in the velocity u(x,t). 
The same result is found and compares well with simulations in the HS case. 


5 Spatial correlations and their effects in the HCS 


Assuming space translation invariance, which is certainly valid in the HCS, we can write 
a hierarchy of equations for the spatial correlations defined as Ck, P = (vj tP Vj + k )P ) at a 
distance of k sites, at time p: 


Co,p+i 

— Cb,p + (cF — 1 )L 1 (Co,p — Ci )P ), 

(11) 

C'ljP+i 

= Ci iP + (1 — a 2 )L 1 (C'o, p — C\ jP ) + (1 + a)L 1 (C , 2 lP — C\ iP ), 

— Ck, p + (l + a)L (Cfc_|_i !P + Ck~i )P 2 Ck,p)i 2 <k< ^ , 

(12) 

C/c,p+l 

(13) 

Cl± i p 

2 

= c-^i 

2 

(14) 


A striking consequence of momentum conservation is the sum rule Co iP +2 ^k, P = 

0. Then we expect that correlations are of the order 0(L -1 ). For example, in the elastic 
limit a = 1, their equilibrium value is (vjVj+i) = —T{L — 1) _1 , V/ ^ 0. We take (11) 
and (13) in the continuum limit, x = (k — 1)/L and t = p/L 3 , and retain only terms up 
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Figure 2: Rescaled velocity profile maximum U ma x = U{x\i,t) as a function of time, where 
Xm = 1/4. Trajectories start from a sinusoidal average velocity profile rt(:r,0) = rtosin(27rx) 
(here Uq = 0.1), which gives hydrodynamic predictions U(x,t) = Uq sm(2'Kx)e^~ l ' 0 ' >t ^ 2 (drawn as 
solid lines). 


to 0(L *), obtaining 


d t D(x,t ) = 2 d xx D(x,t) + 0{L~ 2 ). 


(15) 

(16) 


Here, we have introduced the notations D(x,t) = LC(x,t) and t/(t) = lim x _K) D(x, t). 
Expression (15) introduces a correction in the hydrodynamic average granular tempera¬ 
ture, given by the nearest-neighbour-particle velocity correlation, whereas (16) is a dif¬ 
fusion equation for spatial correlations. Boundary conditions stem from (12) and (14), 
which give a reflecting boundary at x = 1/2 and the sum rule for momentum conser- 
vation, T(t) + 2 f 0 D(x, t ) = 0. In the long time limit, we obtain the following scaled 
stationary solution 


D(x) 


—Hcos 


tt\/-{1-2x) 



(17) 


where D(x) = lim^oo D(x,t)/Tijcs(t). Note that the Fourier transform of D(x) + e)(xrj 
takes the form S(k) = with & = 2vrn and n is a positive integer, which has been 

derived for the correlations of the velocity shear mode in d > 1 from Landau-like granular 


: The delta function is needed to include the case of the autocorrelation (vf), since D(0) corresponds 
to {ViVi+i), see the paragraph above (151. 
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fluctuating hydrodynamics, see for instance 26 . This result reinforces the motivation 


of our Id model as a simple picture for the shear mode of the velocities in d > 1. 

In Fig. [3j we compare the theoretical prediction in ( |17[ ) for the MM case with nu¬ 
merical results. Remarkably, such prediction compares well also with HS simulations, 
where the analytical computation appears to be more challenging. In conclusion, the 
mechanism that induces spatial correlations in the system seems to be independent of 
the particular interaction model. 



Figure 3: Amplitude A of the rescaled correlation function defined in (17) as a function 
of the dissipation rate u, both for Maxwell molecules (MM) and hard spheres (HS). 
Trajectories start from a homogeneous mesoscopic velocity profile u(x, 0) = 0. The 


theoretical prediction of (17) is also shown (line). 


Equation (15) suggests that the Haff law has a finite size correction. We consider a 
perturbation around the HCS T(t) = Tncsit) + L~ 1 6T(t ) + ... and D(x,t ) = Dhcs + 
L~ 1 5D(x,t ) + .... Making use of (15) and defining 5T(t ) = 5T{t)/T hcs{ t), we obtain 
AST = v^hcs, with 


IpHCS 


ipHCsjt) 
THCs(t ) 




Hence, the granular temperature follows 


T(t) = T HCS {t) 


1 + Z ^HCSVt + 0(L 2 ) 


(18) 


(19) 


which is valid for not very long times (t not scaling with L). There is a critical dissipation 
value = I'c/d = 2vr 2 where ipHCS changes sign, and this determines a change of the 
time-derivative of 5T. Thus, at finite (large) values of L , the temperature decays faster 
or slower than the Haff law if v < up, or v > respectively. In Fig. [4j we compare 
the predicted Haff law finite size effect with the simulation results, obtaining excellent 
agreement. 
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Figure 4: Rescaled temperature ST as a function of time, for v = 10, 20,, 70 (from bottom to 
top). Numerical values (circles) are plotted together with linear fits (lines) made on the second 
half of the trajectory. In the inset, we show the comparison between the fitted slopes in the main 
panel (points) and the theoretical ones from (19) (blue line), as a function of v. The horizontal 
black line at i^hcs = 0 shows the transition at i /0 = 2i r 2 . 


6 Summary. 


In conclusion, we have discussed the rigorous hydrodynamic limit of a lattice model for 
granular and active fluids with momentum conservation and energy dissipation. Macro¬ 
scopic equations reproduce the realistic evolution of the velocity shear mode, which is 
diffusive, as well as that of the temperature field, which includes heat diffusion, inelastic 
dissipation and viscous heating. A crucial phenomenon of inelastic fluids, that is, the 
shear instability of the homogeneous cooling state, is recovered. 

The model allows us to derive the evolution of ~ 1/L spatial correlations, which 
present non-trivial long range extension due to momentum conservation and alter the 
temperature decay in an observable way. This opens new interesting paths of investi¬ 
gation, such as trying to relate the deviation from the Haff law found here with the 
renormalization of the cooling rate found in systems of smooth HS near the shear insta¬ 
bility 23 . The appearance of long-range correlations under non-equilibrium conditions 
for conserved fields is a feature well expected on general grounds [5,6], but rarely derived 


in full analytical detail. 

Finally, we stress the importance of considering finite size effects in granular systems, 
since large sizes are rarely realized in experiments. In addition, the HCS is unstable for 
large L , and therefore finite size corrections cannot be disregarded by considering an 
arbitrarily large system. These facts, together with the scarcity of studies about finite 
size effects in granular matter, makes worthwhile deepening this point in the next future. 
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A Details for the derivation of the hydrodynamic equations 

Average fields (over initial conditions and noise realizations) are defined as ui tP = (ty p ), 
Ei tP = {vf p ) and T^ p = Ej p — uf p . The microscopic equations for the evolution of 
averages at time p at site l are obtained by averaging equations ([3]) and Q of the main 
text, obtaining 


u i,p+ 1 u i,p — ( 3i,p ) T ( ji—i,p)i 

Ei t p -\.i — {Ji,p} T (Jl-i,p) T {di t p) • 


( 20 ) 

( 21 ) 


Averages of currents and dissipation can be computed assuming the Local Equilibrium 
approximation (LEA), which is explicitly stated as 


P2(vi,v l+ i;p) 


\/2ltT l) pyj2TtTl + l . j 


( v i~ u i,p) 

2T i,p 


2T l+l ,p 


( 22 ) 


In the Maxwell molecules case (/3 = 0), where one has (S y i) = 1/L, computations using 
the LEA givcj^] 


( di,p) 


(Ji,p) 


( Jl,P ) ~ 2L ( U l,P U l+ l,p) > 

(Ti,p ~ Ti+ i,p + u ^ p ~ u ‘i+i,p ) > 


(23a) 

(23b) 


a 2 — 1 

2 L 


Ti, P + 


T+i,p+T-i,p 

2 


+ 



V;-j-i J p+'U;_i i p 

2 


2 

+ 


'U'l — l,p 

V 2 


2 


.(23c) 


The hydrodynamic limit consists in a change of spatial and time scales, from (l,p) 
to (x,t), related by size-dependent factors: 

x = l/L, t = p/L i . (24) 

This implies for a generic function fi iP 

fl+i,p ~ fl,p = ® (i?) > ( 25 ) 

Ap+i - fi, P = ^r^/OM) + 0 (zff) > (26) 

2 Note that, in the MM case, for obtaining the averages in ( |23| ) the LEA is only used to write that 
(vi,pVi±i,p) = Ui, p Ui± i iP , that is, we assume that velocities at adjacent sites are uncorrelated. This 
hypothesis is somehow similar to the molecular chaos assumption when writing the Boltzmann equation 
for a low density fluid. 
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which are introduced in (20) and (21) to get the final continuous equations in (x,t). 
Each discrete spatial derivative introduces a L _1 factor in the HL. Then, the difference 
between the current terms in the balance equations is of the order of L~ 3 , because the 
average currents {ji. p ) and (Jj iP ) are of the order of L~ 2 . Those terms, therefore, perfectly 
balance the 1/L 3 dominant scaling on the left-hand side, i.e. the time-derivative. Since 
(di }P ) is of the order of (1 — a 2 )/L, to match the scaling 1/L 3 of the other terms, we 
define the cooling rate to be v = (1 — a 2 )L 2 , which is assumed to be order 1 when the 
limit is taken. This choice automatically implies that when L increases one has that a 
approaches unity, a further reason to expect the validity of the LEA. 

By retaining only the highest order terms in the equations, we get expression 0 of 
the main text. It is interesting to note that our expansion in terms of L _1 is similar in 
spirit to the Chapman-Enskog expansion up to Navier-Stokes order, since we are keeping 
up to terms of the second order in the gradients (of the order of k 2 , being k the wave 
vector, in Fourier space). From a purely mathematical point of view, ([T]) becomes exact 
in the limit L —> oo, but v = (1 — a 2 )L 2 of the order of unity, as stated in the previous 
paragraph. Interestingly, the dissipation field d^ p in (23c) admits an expansion in even 
powers of the gradients, as is also the case of granular fluids 37,38 . However, in the 


above limit, the first terms in d; iP including the gradients are of the order of L~ z as 
compared to the contribution —uT at the Navier-Stokes order, that is, they would only 
be considered at the so-called Burnett order. 

From a physical standpoint, 0 is approximately valid whenever the terms neglected 
upon writing it are negligible against the ones we have kept. Since the correlations 
(viVi± i) are expected to be of the order of L -1 as compared to the granular temperature^] 
we must impose that L> 1 and also t <C L. On the other hand, the term proportional to 
the correlations in the evolution equation for the granular temperature is therefore of the 
order of (1—a 2 )L _1 , which must be negligible against the second spatial derivative terms, 
of the order of L~ 2 . Then, (1 — a 2 )L <C 1 must be further imposed when the correlations 
are neglected in 0. This condition, although less restrictive that 1 — a 2 = 0(L 2 ), also 
reinforces the validity of the LEA. On the other hand, when the correlations are fully 
taken into account, as is the case of equations (15) and (16) of the main paper, the value 
of a is not restricted since the only assumption for writing them is that of homogeneity. 


B Computation of the correlations of the hydrodynamic noise 

Noises with respect to averages appear in the currents ji iP = ji :P + £; ;P , J/ )P = J; iP + r)^ p , 
and in the dissipation d/ iP = d; )P + 0; iP , with noises £; iP , r/; iP and 0/ iP defined to have zero 
average. The idea is that each term x is made of a x part which is an average over the 
fast noise (that is, the collisions, which are counted by the fast stochastic variable yi )P ), 
but at fixed vi )P whose evolution is assumed to be slower than noise. 

To obtain the correlations of noise, we exploit a series of conditions. Explicit 
calculations are discussed here for the case of the momentum current noise £/ jP . It 

Tor example, in the elastic case, the correlations {viVi+k) do not depend on the distance k in equi¬ 
librium, and therefore ( ViVi+k) = —T(L — l) -1 , Vfc ^ 0. See Section [H] for more details. 
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is clear that the definition ji p = (1 + a)Ai p /2L corresponds to the above prescrip¬ 
tion for the noise. First, it is straightforward that (6,p6',p') = 0 for p / p' , be¬ 
cause y p and y p i are independent random numbers. Second, we take into account that 
(Sy Pl iSy p> i>) = Sij>(8y pi i) = 5 iji/L, and the fact that all the other contributions are of the 

order of L~ 2 . Thus, for p = p' we have (£/,p£p,p) = (1 + a) 2 S\y /AL + 0(L~ 2 ). 

At this point, the quasi-elasticity of the microscopic dynamics makes it possible to (i) 
substitute (1 + ot)/2 by 1 and (ii) calculate (A f p ) by using the LEA, to obtain 


(6,p&v) ~ £ 2 Ti t pdi t ii 5 p y. 


(27) 


In the large size system, j^ p scales as L 2 (see 32 ). Therefore, the mesoscopic noise of 


the momentum current is defined as £(x, t) = lirri r.^oc T 2 £z iP , and j (x,t) = j(x,t)+£(x,t), 
in which, again, j(x,t ) = lim^oo L 2 ji tP . Going to the continuous limit, remember¬ 
ing (24), and taking in account that Sijr/Ax ~ 5(x — x 1 ) and 5 p y / At ~ 5 (t — t') we get 


the noise amplitude of the velocity current in the main text. Identical considerations 
lead to the amplitude for the energy current noise. For the fluctuations of dissipation, 
the dissipation term is split again as di p = d^ p + with (di }P ) = (d^ p ). We know 
from the dissipation current definition that (Oi^Oi'y) = 0 for p = p'. Making use of the 
LEA and in the large size system di )P scales as L~ 3 and it is expected that the noise 


should have the same scaling. Going to the continuous limit and taking in account (24) 


and (27), the result in the main text is recovered. 


The cross correlations between different noises are straigthforwardly obtained, along 
similar lines, with the result 


(£(x,t)rj(x,t)) = (r](x,t)€(x,t)) = 4T ( X ’*) , 

i€{x,t)9(x,t)) = {0(x,t)£(x,t)} = 0, 

(??(x, t)9(x, t)) = (0(x, t)r](x, t)) = 0. (28) 


C Numerical comparison for the amplitude of hydrodynamic noises 

A comparison for the amplitudes of noise for the velocity and energy currents is shown 
in Fig. [5j A case with the MM interaction (/? = 0) is considered. The simulations are 
performed with periodic boundary conditions, therefore without energy injection, and 
starting with a non-homogeneous initial condition. The initial mesoscopic velocity profile 
and homogeneous granular temperature are u(x,0) = uq sin(27rx) and T(x,t) = To, 
respectively, with uq = Tq = 1 . 
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